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Abstract 

The influence of the isospin- independent, isospin- and momentum-dependent equation of state 
(EoS), as well as the Coulomb interaction on the pion production in intermediate energy heavy ion 
collisions (HICs) is studied for both isospin-symmetric and neutron-rich systems. The Coulomb 
interaction plays an important role in the reaction dynamics, and strongly influences the rapidity 
and transverse momentum distributions of charged pions. It even leads to the tt~ /ir + ratio deviating 
slightly from unity for isospin- symmetric systems. The Coulomb interaction between mesons and 
baryons is also crucial for reproducing the proper pion flow since it changes the behavior of the 
directed and the elliptic flow components of pions visibly. 

The EoS can be better investigated in neutron-rich system if multiple probes are measured 
simultaneously. For example, the rapidity and the transverse momentum distributions of the 
charged pions, the ir~ /ir + ratio, the various pion flow components, as well as the difference of 
7r + -7r~ flows. A new sensitive observable is proposed to probe the symmetry potential energy 
at high densities, namely the transverse momentum distribution of the elliptic flow difference 

[At£ + -*"(pr-)]- 

PACS numbers: 24.10.Lx, 25.75.Dw, 25.75.-q 
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I. INTRODUCTION 



Intermediate energy heavy ion collisions (HICs) are closely connected to the investigation 
of the nuclear equation of state (EoS). One of the main issues is to pin down the incom- 
pressibility (K NM ) of nuclear matter. Although many efforts have been pursued in the past 
:'ew decades this problem is still far from beingsolved thoroughly (see, for example, Refs. 

QQ). 

furthermore, the symmetry energy 

Blag 

::or isospin-asymmetric matter, the 
momentum dependence of the EoS, have also been found to be very im- 

portant for the dynamics of the intermediate energy HICs and make the problem of probing 
the EoS even more complex. 

In order to explore explicitly the incompressibility, the isospin dependence and the mo- 
mentum dependence of the EoS, sensitive probes were put forward, however, mostly indi- 
vidually. Nevertheless, sensitive probes are not always proprietary, that is, they might be 
affected not only by single physical quantity. Therefore, it is quite necessary to explore 
multiple probes simultaneously so that the comparison between the experimental data and 
the corresponding theoretical predictions becomes more consistent. 

In addition to free nucleons and light fragments, probes related to pions such as the 
tt~/tt + ratio and the directed and elliptic flows, have be en p roven to be very useful to 



r s, nave De en p roven to De very i 

HQ 0, Q Q 0, Q Q 0, 



test the reaction dynamics as well as the EoS [1( ilia ua, im 113, Ha im na Ha i2oi, i2i|. 

It was supposed and observed in Refs. 0, IliA I22I] that pions show a weak positive flow 
effect in central collisions while a weak antiflow effect in peripheral collisions due to the 
shadowing effect of spectators. It should be mentioned that experimental measurements on 
the pion production with various collision systems and beam energies have been released by 



the FOPI/GSI Collaboration in recent years, for example, see Refs. 23|,|2J|, which indeed 
largely deepen our insight into the production mechanism of pions in intermediate energy 
HICs, as well as the dynamics of the nuclear reaction itself. 

In this work we will study a variety of observables related to pion production in isospin- 
symmetric and -asymmetric systems from central and semi-peripheral intermediate-energy 
HICs. The microscopic transport model - ultrarelativistic quantum molecular dynamics 



(UrQMD) |25|, l26l 1271 1281 ] - is adopted with an update of the potentials in the mean field 
part. 

The paper is arranged as follows. In section II, we briefly introduce the UrQMD model 
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and the improvements of the potentials in the mean field part. In section III, the results of 
the ratios between the yields of charged pions, the various directed and elliptic flows of pion 
mesons are shown and discussed. Finally, a summary is given in section IV. 



II. THE URQMD TRANSPORT MODEL AND THE POTENTIAL UPDATE 



301 ] . It was shown 



In the UrQMD transport model [2Jj l2f|, the initialization of projectile and target nuclei, 
the equation of motion of hadrons, and the collision term are described microscopically. This 
model was designed at the beginning for simulating HICs in the energy range from SIS to 
RHIC, where the contribution of nuclear mean field potentials to the dynamics of the reaction 
is considered to be weak. It is known that the conventional (or isosp in-dependent) quantum 

, |29( is mainly applied to the (isospin-asymmetric) 
intermediate-energy HICs. The UrQMD model inherits the basic treatment of the baryonic 
equation of motion in the QMD model, thus, after introducing some modern ingredients 
for the mean field part, it is believed that the UrQMD model can also be used to properly 
describe the physical phenomena in HICs at intermediate energies. 

Furthermore, it has been seen that the contribution of the mean field to the reaction 
d—can not sim p. y be negated at SIS and even AGS ene rgl e S Q 
in Ref. [30] that the experimental elliptic flow as a function of beam energy ranging from 
2 — 10 A GeV for midrapidity protons in Au+Au collisions can be only successfully described 
if mean field potentials are taken into account, while pure cascade calculations fail. In Ref. 
Q| it has been shown that the potential that £ hyperons encounter in the nuclear medium 
plays an important role in their evolution process in HICs at SIS energies. Thus, we have 
updated the potential interactions in the UrQMD model for the intended studies. 

In addition to the isospin- and momentum- independent terms originally implemented in 
the mean field part of the UrQMD model, the following contributions are also considered 
here: 

1. The contribution of the Coulomb interaction between mesons and hadrons (mesonic 
Coulomb potential) are supplied in addition to the Coulomb interaction between two baryons 
(baryonic Coulomb potential), in total, we call it hadronic Coulomb potential. It has been 
found that the Coulomb interaction between baryons and charged 7r mesons plays an essential 
role in the dependence of the n~ /n + ratios on rapidity and transverse momentum \v\. I23I . 
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FIG. 1: Parametrizations of the nuclear symmetry potential energy DDHp*, DDH3p5*, F15, and 
Fa3, as well as a linear one as a function of the reduced density u. 

2. Like in our previous work 17, Q], the symmetry potentials of all baryons, i.e., the 
nucleons, the A(1232)s, the iV*(1440)s, and the hyperons A and E, are introduced. Four 
density-dependent parameterizations for symmetry potential energy are considered: (1) u 1 
with 7 = 1.5 (called F15). Here u = p/po is the reduced nuclear density; (2) w(a — u)/ (a — 1) 
with a = 3 (Fa3). a is the so-called reduced critical density 

H; (3) and (4) so-called DDHp* 
and DDH3p5* symmetry potential energies, which are inspired by the relativistic mean-field 



calculations of DDHp and DDH3p5 



321 ] . The symmetry energy coefficient So = 34 MeV 



is adopted |1|, |33[ . The density dependence of the symmetry potential energies adopted is 
shown in Fig. ^ One sees from Fig. [I] that for reduced densities u < 1, DDH3p<5* is very 
close to Fa3, both lie between DDHp* and F15. For u > 1, their density dependences are 
rather different: for 1 < u < 2.6, the order of the symmetry potential energy E£°^ is F15 > 
DDH3p5* > Fa3 > DDHp*, while for u > 2.6, F15 > DDH3p5* > DDHp* > Fa3. 
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TABLE I: Parameter sets for the nuclear equation of state used in the extended version of the 
UrQMD model. 



EoS 


a [MeV] p [MeV] 


7sky 


t md [MeV] 




So [MeV] 


H 


-165 


126 


1.676 






34 


S 


-353 


304 


7/6 






34 


HM 


-138 


60 


2.08 


1.57 


500 


34 


SM 


-393 


320 


1.14 


1.57 


500 


34 



3. Momentum-dependent interactions for all baryons are introduced. The form of the 
momentum dependence is taken from the IQMD model hJQ|> which reads 



U md = ^mdln 2 [l + a md (Ap) 2 ]«, 



(1) 



in which Ap = Pi — Pj represents the relative momentum of two nucleons i and j. The 
parametrizations of t m d and listed in Tabled We note that B.-A. Li et al. have used 

an isospin-dependent momentum-dependent parametrization in the BUU model, which is 
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361 ] . However, 



guided by a Hartree-Fock calculation using the Gogny effective interaction 
here we do not consider the isospin dependence in the momentum dependent part of the 
mean field. 

Besides the two- and three-body Skyrme potential (t/gky = au + /5M 7sk y), similar to the 
IQMD model we consider a variant nuclear incompressibility, namely, a hard EoS 

with K = 300 MeV ("H-EoS"), a soft EoS with K = 200 MeV ("S-EoS"), a hard EoS 
(K = 380 MeV) with momentum dependence ("HM-EoS"), and a soft EoS (K = 200 MeV) 
with momentum dependence ("SM-EoS"). The parameters of the various EoS used in the 
UrQMD model are listed in Table |U In order to reproduce the ground-state properties 
of finite nuclei in the UrQMD model (for instance, the binding energy Eb and the r.m.s 
radius), the parameters have been slightly readjusted. Note, that the potentials used here 
are not included in the currently available version of the model, but will be included in 
future versions. 
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III. RESULTS AND DISCUSSIONS 



Firstly, we test the applicability of the UrQMD model to intermediate energy HICs. Fig. 
121 shows the results for the excitation function of pion multiplicities in central (6 = fm) 
197 Au+ 197 Au reactions at energies from OA A to 1.5 A GeV and the comparison with the 



recent FOPI preliminarily experimental data 



3q . The ratios between the calculations (with 



S- and SM-EoS and without potentials, namely the cascade mode) and the experimental data 
are shown in the lower plot. In Ref. |39| the multiplicities of pions were compared within 
various transport approaches for HICs around \A GeV. A significant variation of the total 
pion yield with the different approaches has been found. Generally, these transport models 
overpredict the pion yields. Overall, our calculations (Fig. |2] with EoS, especially with a 
S-EoS) are in good agreement with data, while the results without mean-field potential (the 
cascade mode) deviate from the experimental data at lower beam energies (< 0.8A GeV). 

Now let us turn to see how both the Coulomb and the symmetry potentials influence the 
production of pions in isospin-symmetric HICs. Fig. El shows the rapidity (upper plots) and 
the transverse momentum (lower plots) distributions (in the center-of-mass system) of pions 
in central (6 = fm) collisions 40 Ca+ 40 Ca at beam energy Ef, = OAA GeV. The SM-EoS and 
the F15 symmetry potential energy are selected here. In the left (right) two plots we show 
the distributions without (with) the contribution of both Coulomb (of all charged hadrons) 
and symmetry potentials. When the Coulomb and symmetry potentials are switched off, 
the multiplicities of 7r~, 7r + , and 7r° are equal to each other within statistical errors, which 
means that the isospin symmetry is preserved, as one expects. While the Coulomb and 
symmetry potentials are switched on, the rapidity and transverse momentum distributions 
of 7r~, 7T + , and 7T° mesons are clearly different. At mid-rapidity and in the low transverse 
momentum region, the multiplicity of tt + mesons is smaller than the ir~ ones, and vice versa 
in the region of the projectile-target rapidity or large transverse momentum. We notice 
that, after integrating the multiplicities of n~ and 7r + mesons separately, the n~ /ir + ratio is 
about 1.05. We have also checked this ratio at a higher beam energy, = 0.8A GeV, and 
found that it is reduced to 1.02. 

Fig. |U shows the rapidity (upper plot) and the transverse momentum (lower plot) dis- 
tributions of the 7r _ /7r + ratio for the same reaction as in Fig. El In the plots we show the 
results for the following cases (1) without hadronic Coulomb and baryonic symmetry po- 



7 



100 



5^ 

o 

'-4— ' 

E 

a 
o 
'o. 



10- 



Ai i+Ai i fa) 


h=0 fm 

u — \J 1 1 1 1 












5 




I * 


5 










t 






4 
















* 




■ 


FOPI Data 


• 




• 


SM-EoS 


T 




A 


S-EoS 






T 


Cascade 



0.4 



0.6 



0.8 1.0 
E h [A GeV] 



1.2 



1.4 



1.6 



2.0 
1.8 
1.6 
1.4 

- 1.<M 

UJ 0.8H 
0.6 
0.4 H 
0.2 
0.0 



Au+Au (q) b=0 fm 




0.4 



0.6 



Exp/SM-EoS 



- - R 

..a.. R 



Exp/Cascade 
Exp/S-EoS 



0.8 



1.0 



1.2 



1.4 



E b [AGeV] 



1.6 



FIG. 2: Upper plot: The excitation function of the pion multiplicities for central (6 = fm) 
Au+ 197 Au reaction. The Cascade- mode, S-EoS, and SM-EoS are adopted for calculations. The 
FOPI preliminary data for central Au+Au collisions are also shown 38]. Lower plot: The ratios 
between the experimental data and the calculations. 
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tentials ("No Coul. and Sym. Pot.") (2) with hadronic Coulomb and baryonic symmetry 
potentials F15 and DDHp* ("F15" and "DDHp*") and (3) without the hadronic Coulomb 
potential but with the DDHp* symmetry potential energy ("No Coul. Pot."). One finds 
that without the Coulomb and symmetry potentials the it~ /tt + ratio is around unity. When 
the Coulomb and symmetry potentials are switched on, the it~/tt + ratio depends weakly 
on the symmetry potential. At mid-rapidity and at large transverse momenta, the softer 
the symmetry potential is, the larger the tt~/tt + ratio is, and the other way around in the 
regions of projectile-target rapidity and small transverse momentum. When the Coulomb 
interaction is switched off but the symmetry potential is on, i.e., case (3), the tt~ /tt + ratio 
is again around unity. This means that the Coulomb interaction is the leading cause of the 
7r~/7r + ratio deviating from unity, which is shown in case (2). When the Coulomb potential 
is taken into account, the protons and also the positively charged As are slightly pushed 
into the lower density region so that more neutrons and negatively charged As will be in 
the high-density region. As a consequence, the isospin asymmetry (p n — p p )/(p n + Pp) will 
be different from zero locally although the total system is isospin-symmetric. Thus, the 
symmetry potential begins to play a role that leads to the 7t~/tt + ratio depending on the 
symmetry potential for case (2). So we expect that the deviation of the n~ /ir + ratio from 
unity depends on rapidity and transverse momentum and helps to determine the symmetry 
potential for isospin-symmetric systems. 

Before investigating the influence of various mean-field potentials on the pion flow in 
intermediate energy HICs, it is necessary to show the rapidity distributions of 7r~ and 7r + 
mesons (Fig. EJ), as well as the tt~/tt + ratios (Fig. IHJ) ■ They are calculated (1) with different 
isospin-independent EoS, that is, the H-, S-, and SM-EoS (left plot, titled "EoS-0") but 
with the same symmetry potential F15 and including the hadronic Coulomb interaction; 2) 
with or without mesonic Coulomb potential (middle plot, titled " M-Coul-Pot" ) but with the 
same H-EoS and the symmetry potential F15; and (3) with different symmetry potentials 
F15 and Fa3 (right plot, titled "Sym-Pot") but with the same H-EoS and the hadronic 
Coulomb potential. Semi-peripheral (6 = 7 — 9 fm) 208 Pb + 208 Pb collisions at Eb = 0.8A 
GeV are studied for all following calculations. From Fig. |S] we can see, at first glance, the 
uncertainty of the isospin-independent EoS obviously affects the total pion multiplicity; the 
effect of the density dependence of the symmetry potential on the 7r~ multiplicity is visible, 
while the mesonic Coulomb potential hardly affects the multiplicity of charged pions in the 
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FIG. 3: The rapidity (upper plots) and the transverse momentum (lower plots) distributions of 
pions with (right plots) or without (left plots) hadronic Coulomb and baryonic symmetry potentials. 
The isospin-symmetric reaction 40 Ca + 40 Ca at the beam energy Ej, = OAA GeV with the impact 
parameter b = fm is chosen. 
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FIG. 4: The rapidity (upper plot) and the transverse momentum (lower plot) distributions of 
the 7r _ /-7r + ratio as calculated (1) without symmetry and the Coulomb potentials ("No Coul. and 
Sym. Pot."); (2) with symmetry and Coulomb potentials (cases "F15" and "DDHp*"); and (3) 
without Coulomb potential but with symmetry potential DDHp* ("No Coul. Pot."). The reaction 
40 Ca + 40 Ca at E b = OAA GeV with b = fm is chosen. 
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FIG. 5: Rapidity distributions of it and tt + for various isospin-independent EoS (S-, H, and SM- 
EoS, left plot), with or without mesonic Coulomb potential (middle plot), and with the F15 and 
Fa3 symmetry potentials (right plot) (see text). The reaction 208 Pb + 208 Pb at E b = 0.8A GeV 
and b = 7 — 9 fm is chosen. 

neutron-rich 208 Pb + 208 Pb system. These characteristics have already been seen in jl7|. 
The left plot of Fig. El shows that the pion yield with momentum dependent interactions 

SM-EoS) is obviously lower than that with S-EoS or H-EoS, which has been shown in Ref. 

3|. Further, we see that both the mesonic Coulomb and the symmetry potentials hardly 
affect the 7r + multiplicities in neutron-rich HICs. Therefore, the rapidity distribution of 7r + 
can provide more accurate information on the isospin-independent part of the EoS for the 
neutron-rich HICs. 

The 7r~/7r + ratio at — 1 < y c < 1 is smaller for a soft EoS as compared to a hard EoS (left 
plot of Fig. EI), due to more two-body scatterings. The 7r~ /tt + ratio with a soft symmetry 
potential is larger than the one with a hard symmetry pot ential (right plot). The effect of 
the symmetry potential has been studied extensively in [itL ll^ . Here we want to stress that 
both, the isospin-independent and isospin-dependent EoS parts, obviously affect the n~ /n + 
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ratio but in the opposite direction when the stiffness changes from soft to hard. When the 
mesonic Coulomb potential is not taken into account (middle plot) the tt~/tt + ratio is larger 
at rapidities exceeding one than at mid-rapidity, which was proven wrong by experimental 
measurements, for example, see Ref. 24]. It means that the mesonic Coulomb potential 
needs to be considered, besides the baryonic Coulomb potential. One further finds from this 
figure that the tt~/tt + ratio becomes rather flat at small rapidities (—0.75 < y c < 0.75) in 
the semi-peripheral HICs. We have also found a double peaked structure (in the projectile 
and target regions) of the rapidity distribution of the ratio at even larger impact 

parameters. This behavior is rather different from central collisions (upper plot of Fig. EJ) 
where the ir~ /tt + ratio is peaked in the mid-rapidity region. This change of the behavior of 
the 7r~/7r + ratio as a function of rapidity from central to semi-peripheral HICs is due to the 
decrease of the number of nucleon-nucleon collisions 

Now let us turn to discuss collective flow observables. Fig. [7| shows the in-plane directed 
transverse flow of charged pions as a function of rapidity. It is well known that the directed 
flow of particles is particularly sensitive to the stiffness of the EoS, which is also shown in 
the left plot of Fig. The weak anti-flow of charged pions for semi-peripheral collisions 
is well reproduced, similar to the calculations and the experiments in Refs. 0, Q 
for semi-peripheral Au+Au collisions at different beam energies. Among the three types of 
the EoS studied here, the anti-flow effect of ix~ and tt + calculated with the SM-Eos is the 
strongest, followed by the S-EoS jlOj. The anti-flow effect of tt + is much stronger than that 
of 7r~. If the mesonic Coulomb potential is neglected, the absolute value of the antiflow 
parameter of tt~ is even higher than that of 7r + as shown in the middle of Fig. which is in 



disagreement with experimental measurements [22j . In the right plot of Fig. [71 we find that 
the effect of the density dependence of the symmetry potential on the charged pion flow is 
very weak, which means that the in-plane directed flow of charged pions can provide more 
accurate information on the isospin-independent part of the EoS rather than on the isospin 
dependent part. 

In Fig. Elwe show the rapidity distributions of the 7r + -7r~ transverse flow difference Ap^. m 
defined by Ap^ m = p 1 ^ — p^ . p^ m is similar when calculated with S-EoS, H-EoS as well 
as SM-EoS at mid-rapidity, but it turns out to be different between with the momentum 
dependence (SM-EoS) and without the momentum dependence (S-EoS and H-EoS) in the 
projectile-target rapidity region. Again, from the right plot we find that ApV m depends 
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FIG. 6: Rapidity distributions of it~ /tt + ratios for various isospin-independent EoS (S-, H, and 
SM-EoS, left plot), with or without mesonic Coulomb potential (middle plot), and with the F15 
and Fa3 symmetry potentials (right plot) (see text). The reaction 208 Pb + 208 Pb at Ej, = 0.8A 
GeV and b = 7 — 9 fm is chosen. 

only weakly on the symmetry potential. Thus, we conclude that the 7r + -7r _ transverse flow 
difference Ap^. m can be adopted to probe the momentum dependence of the EoS. 

Fig.Elshows the influence of different parts in the EoS on the ir + -ir~ elliptic flow difference 
as a function of rapidity, which is defined as Auf™ = v% + — fj • Here t>2 = ((ff) 2 — (^f ) 2 ) anc ^ 
t>2 + is the elliptic flow for 7r + , v J is the elliptic flow for n~. First, we see that Av^™ value is 
always negative when the mesonic Coulomb interaction is taken into account. The influence 
of the mesonic Coulomb interaction on the 7r + -7r~ elliptic flow difference is shown in the 
middle plot of Fig.^l Av 2™ is almost zero when the mesonic Coulomb interaction is switched 
off. When the mesonic Coulomb interaction is switched on, the out-of-plane elliptic flow 
becomes larger for 7r + than for 7r~ mesons which leads to a negatively elliptic flow difference. 
The left and right plots of Fig. El show that Av ^ is insensitive to the uncertainties of both 
the isospin-independent and -dependent parts of the EoS. The effect of the momentum 
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FIG. 7: The directed transverse flow distributions of the 7r~ and ir + mesons for various isospin- 
independent EoS (S-, H, and SM-EoS, left plot), with or without mesonic Coulomb potential 
(middle plot), and with the F15 and Fa3 symmetry potentials (right plot) (see text). The reaction 
208p b + 208 pb at Eb = o GeV and b = 7 - 9 fm is chosen. 

dependence of the EoS has only little effect on Av^™- The weak dependence of the rapidity 
distribution of At^" 1 on the mean-field potentials might be useful for extracting exclusive 
information on medium corrections of the binary cross sections. 

From the above discussions on the pion flow (Figs. EKJ), one might have the impression 
that the pion flow is not sensitive to the density dependence of the symmetry potential. 
Actually, this is not true, when we study the transverse momentum dependence of Afg™. 
In Fig. we illustrate Au^" 1 of charged pions as a function of transverse momentum pf m \ 
The SM-EoS is adopted in the calculations. Aff™ is negative for all pf m \ The dependence 
of AfJ? 111 on the form of the density dependence of the symmetry potential becomes very 
pronounced at large transverse momenta ~ 0.2 — 0.5 GeV/c. This behavior is very similar 



to the proton and neutron elliptic flow difference v^ 1 32j,|40j. As it is known that pions with 



high transverse momenta are mainly produced from the high density region, Comparing 
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FIG. 8: The rapidity distributions of the 7r + — 7r~ transverse flow difference (ApP m = p^ —p^ ) for 
various isospin-independent EoS (S-, H, and SM-EoS, left plot), with or without mesonic Coulomb 
potential (middle plot), and with the F15 and Fa3 symmetry potentials (right plot) (see text). The 
reaction 208 Pb + 208 Pb at E b = 0.8A GeV and b = 7 - 9 fm is chosen. 

Fig. [TUl with Fig. [TJ we deduce that it reflects explicitly the density dependence of the 
symmetry potential at densities u > 1. At densities u > 2.6 (Fig.QJ the symmetry potential 
Fa3 is smaller than DDHp* which is also reflected in Fig. El a t the transverse momentum 
~ 0.5 GeV/c. At higher transverse momenta, although more than 1.6 million events were 
calculated, the statistical errors are quite large and the results are not shown here. 



IV. SUMMARY AND OUTLOOK 



In summary, based on the UrQMD model (version 1.3) we have studied the influence of 
different parts of the EoS such as the isospin-independent, the Coulomb interaction (the 
Coulomb interaction between meson and meson (baryon) is also considered) and the sym- 
metry energy on a variety of observables related to pion production in intermediate energy 
HICs. We first have studied the isospin-symmetric central 40 Ca + 40 Ca collisions. It shows 
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FIG. 9: Rapidity distributions of the tt + — ir~ elliptic flow difference (Av^ 111 = v\ — v% ) for 
various isospin-independent EoS (S-, H, and SM-EoS, left plot), with or without mesonic Coulomb 
potential (middle plot), and with the F15 and Fa3 symmetry potentials (right plot) (see text). The 
reaction 208 Pb + 208 Pb at E b = 0.8A GeV and b = 7 - 9 fm is chosen. 

that the Coulomb interaction plays an important role in the reaction dynamics and influ- 
ences the rapidity and transverse momentum distributions of charged pions considerably. 
Due to the effect of the Coulomb potential, the ir~ /tt + ratio deviates from unity and shows 
the dependence on the symmetry potential even for isospin-symmetric systems. Our study 
shows that the rapidity distributions of the pions and the 7[~/tt + ratio are strongly influ- 
enced by the uncertainty of the EoS, while positively charged pions from neutron-rich HICs 
are much less influenced by the density dependence of the symmetry potential. 

We have also studied the pion directed and elliptic flow (difference) as functions of rapidity 
and transverse momentum for the neutron-rich system 208 Pb + 208 Pb. We find that the 
mesonic Coulomb interaction plays an important role in reproducing the proper pion flow. 
The omission of mesonic Coulomb interaction leads to a wrong behavior of the rapidity 
distribution of the n~ /n + ratio and the charged pion directed and elliptic flow. We have also 
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FIG. 10: The transverse momentum distribution of the A^J,™ with different symmetry potentials. 
The SM-EoS is adopted for the 208 Pb + 208 Pb reaction at E b = 0.8A GeV and b = 7 - 9 fm. 

found that the directed flow is sensitive to the isospin-independent EoS but not sensitive 
to the various forms of the density dependence of the symmetry potential. The rapidity 
distribution of the transverse flow difference Apg 111 of charged pions shows the sensitivity 
only to the momentum dependence of the EoS. The rapidity distribution of the elliptic flow 
difference Av^" 1 of charged pions is insensitive to the mean field potentials. However, the 
transverse momentum distribution of Av% m of charged pions (at transverse momenta ~ 0.2 — 
0.5 GeV/c) becomes more sensitive to the form of the density dependence of the symmetry 
potential. Thus, we can use the extensive comparison of multiple observables between 
the calculations and the intermediate-energy experimental data to extract comprehensive 
information on the EoS. 
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